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;h ! Abstract 

' First order quark hadron phase transition is considered in scaling hydrodynamics 

including the Csernai-Kapusta model of nucleation of hadronic bubbles leading to 
supercooling. The effect of viscosity on the entropy production is studied in ideal as 
well as in viscous hydrodynamics with and without supercooling. It is found that 
the excess entropy produced due to supercooling depends weakly on viscosity. 
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^ : 1 Introduction 

o 

' The ultrarelativistic heavy ion colhsions provide a means to create a new state of matter 

at high temperature and density as the quantum chromodynamics (QCD) predicts a 
phase transition from normal hadronic matter to quark gluon plasma (QGP), a state 
of unconfined quarks and gluons. In recent years a considerable amount of theoretical 
and experimental activity is going on in this field. After the collision of two Lorentz 
' contracted nuclei, QGP is assumed to be formed and equilibrated in a very small time 

of the order of ~ 1 fm. In the usual description of the evolution of the plasma, it cools 
^ ! by expanding hydrodynamically till the critical temperature is reached at which a 

^ I transition from QGP to normal hadronic matter takes place. The temperature remains 

fixed at until hadronization gets completed. This hadronization scenario corresponds 
to the ideal Maxwell construction. However, in reality hadronization does not begin at 
T = Tc due to the large nucleation barrier. Recently Csernai and Kapusta have proposed 
a model for nucleation for the relativistic first order phase transition [0]. Supercooling, 
through nucleation of hadronic bubbles in QGP has been studied by several authors 
H, 1^, ^ ^ using the Csernai-Kapusta model of nucleation. A general outcome of these 
studies is that the plasma will cool according to the law T(t) = Tq(tq/tY^^ till T^. The 
matter continues to cool below until the temperature goes down to about ~ .8 T^, where 
bubble formation and growth becomes sufficient to reheat the system due to the release 
of latent heat. Compared to the idealized Maxwell construction the supercooling delays 
the transition and leads to an extra entropy production as the nucleation process allows 
dissipation around the hadronic bubbles. The dynamical prefactor in the nucleation 
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rate includes quark viscosity coefficients which bring dissipative effect in the medium. 
In fact, the nucleation rate is hmited by the ability of the dissipative processes to carry 
latent heat away from the bubbles's surface. However, for the dynamical evolution of the 
plasma the ideal hydrodynamics is used. This is not consistent as the viscosity dependent 
terms in hydrodynamics would also contribute to the entropy production. Therefore, in 
this work we study the supercooling in the viscous hydrodynamics. The role of viscosity 
on supercooling and entropy production has been investigated in detail. 

2 Csernai-Kapusta model of nucleation 

The nucleation model computes the probability that a bubble of the hadronic matter 
appears in a system, initially in QGP phase near the critical temperature. Further in the 
model a small baryon chemical potential is assumed. The bubble formation is assumed to 
take place in a homogeneous QGP phase consisting of u and d quarks and gluons, ignoring 
the role of inhomogenities such as strange quarks. It is further assumed that there is not 
substantial supercooling. Langer's theory of nucleation gives the nucleation rate per unit 
volume at temperature T as 



where AF^, is the change in the free energy of the system due to the formation of a 
critical hadronic droplet, Qq is a statistical prefactor which measures the available phase 
volume and V is the volume of the system. The dynamical prefactor, k determines the 
exponential growth rate of critical droplets which are perturbed from their equilibrium 
radius R^. The coarse-grained effective field theory approximation to QCD is utilised to 
obtain k and flo- The factor is obtained as 



where a is the surface free energy and is determined as 50 MeV/fm^ by lattice guage 
theory simulations without dynamical quarks. The correlation length is estimated as 
0.7 fm. The dynamical prefactor k, is obtained as 



I 



'0 -AF,/T 





(2) 



K — 



(Au;)2i?3 



(3) 



2 



where Acu is the difference in the enthalpy densities of the two phases, rj and ( are 
respectively, the shear and bulk viscosity coefficients. Inserting Qq from Eq. @ and k 
from Eq. in Eq. (|T]) we get the nucleation rate per unit volume as 



-AF,/T (A\ 



The critical radius R^, is given by the Laplace formula as 

R.{T) = (5) 

where Pq/h is the pressure of the quark/hadron phase at temperature T and the AF,, as 

Air 

AF, = -^aRl (6) 

From the nucleation rate Eq. (^), one can calculate the fraction of volume /i(r) which 
has been converted from QCD plasma to hadronic gas at proper time r. If the system 
cools to Tc at time Tc, then at some later time r the fraction h of space which has been 
converted to hadronic gas is 



h{r)= / dr'I{T{r'))[l-h{r')]V{r\r). (7) 

Jtc 

Here V"(r', r) is the volume of a bubble at time r which had been nucleated at an earlier 
time r'; this takes into account the bubble growth. The factor 1 — /i(t') accounts for the 
fact that new bubbles can only be nucleated in the fraction of space not already occupied 
by the hadronic gas. The model for bubble growth is simply taken as p[ 



V{T\r) = ^ (i?.(T(rO) + £ dT"v{TiT")) 
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where v{T) is the velocity of the bubble growth at temperature T. By definition a critical 
size bubble is metastable and will not grow without a perturbation. The growth of bubbles 
has been studied numerically with relativistic hydrodynamics by Miller and Pantano 0. 
Their results are consistent with the growth law 



v{T) = 3c[l - T/T^f'^. (9) 
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This expression is intended to apply only when T > |Tc so that the growth velocity stays 
below the speed of sound of a massless gas, c/-\/3. At the critical temperature, i?* — >• oo, 
AF* oo and the rate of nucleation vanishes. The system must supercool at least ~ 5 
% to attain a finite rate. In the evolution of the matter from QGP to hadron phase, the 
temperature varies with r and its description by scaling hydrodynamics provides another 
equation so that h{T) and T(r) are determined from these equations at any proper time 
r. 



3 Viscous hydrodynamics 

The scaling viscous hydrodynamics is discussed by Danielewicz and Gyulassy 0] and 
others. The form of the dissipative terms depends on the choice of the definition of what 
constitutes the local rest frame of the fluid. The Landau-Lifshitz definition is appropriate 
for describing systems with small (or zero) chemical potential. Here we give a simple 
derivation of the viscous hydrodynamics. In scaling hydrodynamics the expansion takes 
place only along the direction of collision which we chose as z axis. The proper time r 
and space-time rapidity y are used in place of t and z which are defined by 



r = - ^2 (10) 

and 

1, t + z , , 

y=-ln . 11 

^ 2 t-z ^ ' 

In the rest frame of the fluid we take a volume element as 5V = A±T6y, where A± is the 
area of the element transverse to z direction. The expansion takes place with the velocity 
Vz = z/t, in the rest frame of the element. Due to the viscous effects in longitudinal 
direction the heat density per unit time is given by |^ 

or (j) = {Ar]/3 + 0/r^- (12) 

After time increment At, the volume expands by AV = A±AT6y. The amount of work 
done in expansion is pAV. If the energy density at r is e and at r + At is e + Ae then 
the energy conservation implies 

e6V = {e + Ae){6V + AV) + pAV - (pTAV, 
or ~ AeA_iT6y + {e + p- (j)T)A_iAT6y, (13) 
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which leads to the scahng Navier-Stokes equation in the hmit Ar 

^ e + p ^ Ar]/3 + C ^^^^ 

Equation (JT^) has been solved earlier for the QGP to hadron transition with the 
Maxwell construction, i.e., assuming T = for the mixed phase. 

To solve Eq. (|T^ we require knowledge of the equation of state and the temperature 
dependence of rj and (. In this work we use the bag equation of state for QGP. The energy- 
density, pressure and entropy densities in pure QGP and hadron phases are taken as 

eg{T) = 3a,T^ + B, p,{T) = a,T^ - B, s,{T) = Aa^T^ (15) 



eh{T) = 3af,T\ ph{T) = a^T^ s^iT) = Aa^T^ (16) 

Here, and are related to the degrees of freedom operating in two phases and B is 
the bag pressure. 

For ultrarelativistic gases, the bulk viscosity ( is usually much smaller than the shear 
viscosity r] 0. Danielewicz and Gyulassy 0] give the acceptable range of rj for the 
applicability of the Navier-Stokes equation to the expansion of the plasma as 

2T^ <r]< 3T'\tT). (17) 

We define = (4?7/3 + C) and assume the temperature dependence of viscosity coefficient 

as 



with q being replaced by h for the hadronic phase. 

Solving Eq. (p!4D for temperature in pure quark phase, we get 



T = Tn ( ^ 
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where Th is the time where hadronization gets completed. 

The energy density in the mixed phase at a time r can be written in terms of hadronic 
fraction hi^r) as 

6(r) = e,{T) + {e,{T)-e,{T))h{r), 

= 3[a, + {ah - ag)h{T)]T^ + B[l - h{T)], (21) 

while the enthalpy density is 

u{t) = 4[ag + {ah - a,)/i(r)]T^ (22) 

The temperature T can be deduced as 

^^""^ ^ {A[ag+{a!-a,)h{T)]) ' ^^^^ 
The shear and bulk viscosities are also taken to be the functions of time according to 

/i(r) = fiq{T) + (^ft(T) - fig{T))h{T), 
or /i(r) = [figo+{i2ho- l^qo)h{T)]{T{r)f. (24) 

For the Maxwell construction, i.e., for T = T,, in the mixed phase, solution of Eq. (0) 

is 

h{T) = {c- ab)[Et{b/r) - Ei{h/T,)]e-^'^/T + (a + 1)[1 - r./re^^l^^-^'^ (25) 
Here, Ei is exponential integral and 

a = - — - — , 6 = — - — , c= — - — . (26) 

The constants a, b, and c are evaluated for T = Tc. 

The solution given by Eq. (p5|) does not account for supercooling. Due to supercooling, 
i.e., T 7^ Tc in the mixed phase, the temperature is not constant and depends on r. 
Equation (|T^) then becomes 

^ ^ _^ /igo + (/^feo - fiqo)h{^) ^ 

r ^ [4ag + 4{ah-ag)h{T)f/^ ' ^ ^ 
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Equation ( |2^ ) and Eq. (0) are coupled equations, finally to be solved for /i(r) and e(r) 
[ or T(r)]. Once we get /i(r) and T(r) we can calculate the entropy density. 
Eq. ( P7D can also be written in terms of entropy as 

ds ^ s f PhjT) -pg{T) \ dh /Xgo + (/ifep - /igo)/i(r) 

r 1^ T j rft [4a, + 4(aft-a,)/i(r)]2/3 ^2 ' ^ ^ 

Without the second and third terms on the right hand side, this equation describes the 
conservation of entropy. The second term is responsible for the entropy production due to 
nucleation. For the Maxwell construction, p^ — pg = as T = Tc and this term vanishes. 
The third term leads to continuous entropy production due to dissipative effects. We 
discuss the solution of these equations in the next section. 



4 Results and Discussion 

By solving together Eq. (^) and Eq. (0), we have studied the plasma evolution and 
calculated st — the entropy production as a function of r. The initial conditions are taken 
at T = Tc as T = Tc, h = and e(rc) = eg(rc). Now, /i(r) is calculated using Eq. (p5|) 
with step A in r. With this value of h, Eq. (pTf) is solved for e(r) [or T(r)]. Then 
Eq. d^) is evaluated by the trapezoidal rule, thereby yielding new value of h. Using new 
value of h, we solve again Eq. (pTl) to improve the value of e. This is repeated till an 
accuracy of ~ 10~^ is obtained. Then we proceed to the next step Tc + 2A. We take 
aq = 377r2/90 and ah = 4.671^/90. The value 4.6 is taken instead of 3 to account for 
p, uj and f] mesons apart from pions. We chose Tc = 160 MeV in this work implying 
2^1/4 _ MeV. For the hadronic matter r]ho = 1-5, Cho = 1 while for the quark matter 
Cgo = 0, IT^]. The coeficient rjgo has been varied from 2.5 to 20. The nucleation 
involves the surface free energy, correlation length, and the velocity of bubble growth, 
which have already been described in Sec. II. The volume of a fluid element is A±rSy. As 
rSy does not change with r in the scaling hydrodynamics, st provides a measure of total 
entropy. So as to understand the role of supercooling and viscous heat generation, we 
have compared various scenarious; ideal hydrodynamics (IHD), IHD with supercooling, 
viscous hydrodynamics (VHD) and VHD with supercooling. We assumed that at initial 
time Ti = 1 fm/c, the temperature Tj=268 MeV marked the beginning of the evolution. 
This corresponds to energy density e=8.51 Gev/fm^ and entropy Sirj=40.65 fm"^. These 
are appropriate for the Relativistic Heavy Ion Collider (RHIC) energies. For further 
discussions, it is instructive to examine r versus T plot. Figure 1 shows such a plot for 
VHD calculations with supercooling included. From a initial point (rj,Tj), one moves on 
to (tc, Tc) where critical temperature is reached. However, the nucleation rate is zero at 
this point and system does not hadronize. The system keeps cooling to (rm,Tm) where 
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significant nucleation rate and hadronization is reached. The system cannot cool further 
as entropy cannot decrease according to the second law of thermodynamics. At this point, 
the hadron fraction is around 11-18 % for various cases. As the temperature increases 
towards Tc, the hadron fraction h increases and slowly approaches a value of one at the 
point {Tfj_,Tc). After this, temperature starts decreasing finally reaching the point {Tf,Tf) 
at the freeze-out temperature Ty. This is the general character of r — T curve. If one 
uses the Maxwell construction, one reaches directly from (tc, T^) to (t/j, T^) without any 
change in the temperature. 

In Figs. 2 and 3. we show the results for T, h and st as a function of r for riqQ=2.5 
and 5 . For IHD, the critical temperature Tc is reached earlier than VHD. Similarly 
with supercooling included, Tm is approached earlier in IHD than VHD. The supercooling 
leads to an abrupt change in entropy starting from Tm onwards, i.e., in the reheating 
region. We find in our calculations that st does not increase much beyond r ~20 fm/c. 
Therefore, ShTh represents reasonably well the total entropy produced. As Sh is evaluated 
at Tc for all cases, the total entropy is also a measure of the life time of the system upto 
the hadronization stage. Figure 4 shows ShTh as a function of rjqQ. With VHD, entropy 
production increases with viscosity. Supercooling leads almost to a constant shift of the 
curve for the VHD. In Fig. 5, the excess entropy production due to supercooling is shown 
as a function of viscosity. For all values of rjqo, excess entropy does not change much. It is 
almost constant for the VHD calculation. We also found that the values of ri^o << 1 do 
not lead to significant hadronization and the system continues to cool in the QGP phase. 
Sufficient value of rjgo is required for the transition to the hadron phase. The bubble 
growth velocity coefficient (3c in our calculation) also plays an important role in the 
transition as setting it to zero also blocks the transition. The viscosity plays very crucial 
role in the transition through nucleation. However, once viscosity enters into picture, the 
viscous hydrodynamics is to be used for consistency. In the presence of viscosity, the 
initial energy density estimates from the experimental rapidity density distribution would 
be lower than the ones calculated with the Bjorken formula. Though supercooling does 
not lead to significant increase in entropy production, it is very much needed. For the 
description of the phenomena in relativistic nuclear collisions, the viscous hydrodynamics 
with supercooling is a more appropriate framework than the ideal hydrodynamics with or 
without supercooling. Further theoretical work is needed to compute the viscosity with 
narrower bounds for the analysis of the experimental data. 

5 Conclusions 

We have studied the entropy production both in ideal and viscous scaling hydrodynamics 
with and without supercooling. Excess entropy produced due to supercooling in viscous 
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hydrodynamics, weakly depends on viscosity of the plasma phase. Though, in general 
entropy produced due to supercooling is much less than that due to viscous heat gen- 
eration, the phenomenon of supercooling provides an important physical mechanism for 
the quark-hadron transition to occur. The viscous hydrodynamics with supercooling also 
leads to an increase in the lifetime of the plasma. 
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FIG. 1. T/Tc as a function of proper time r in viscous hydrodynamics with super- 
cooling. 

FIG. 2. The T/Tc, h{T) and sr as a function of r in IHD (dotted), IHD with 
supercooling (short dashed), VHD (long dashed) and VHD with supercooling (solid) 
at r]qQ—2.b. 

FIG. 3. Same as Fig. 2, but with 77^0=5.0. 

FIG. 4. ShTh versus 77^0 in IHD, IHD with supercoohng, VHD and VHD with super- 
coohng. 

FIG. 5. Excess entropy produced due to supercooling versus 77^0 both with IHD and 
VHD. 
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